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Abstract 

We present a new determination of the strong coupling constant a s through the 
scaling violations in the fragmentation functions for charged pions, charged kaons, 
and protons. In our fit we include the latest e + e~ annihilation data from CERN 
LEP1 and SLAC SLC on the Z-boson resonance and older, yet very precise data 
from SLAC PEP at center-of-mass energy = 29 GeV. A new world average of 
a s is given. 

PACS numbers: 13.65.+i, 13.85.Ni, 13.87.Fh 



The strong force acting between hadrons is one of the four fundamental forces of 
nature. It is now commonly believed that the strong interactions are correctly described 
by quantum chromodynamics (QCD), the SU(3) gauge field theory which contains colored 
quarks and gluons as elementary particles. The strong coupling constant oiT f (/i) = 
g 2 /(4ii), where g s is the QCD gauge coupling, is a basic parameter of the standard model 
of elementary particle physics; its value (M z ) at the Z-boson mass scale is listed among 
the constants of nature in the Review of Particle Physics JIJ . Here, \i is the renormalization 
scale, and rif is the number of active quark flavors, with mass m q <C fi. The formulation 

of as f \f^) in the modified minimal-subtraction (MS) scheme, with four-loop evolution 
and three-loop matching at the flavor thresholds, is described in Ref. 0. 

There are a number of processes in which ctf\Mz) can be measured (see Refs. |]||, 
for recent reviews). A reliable method to determine af\Mz) is through the extraction 
of the fragmentation functions (FF's) in the annihilation process 

e+e" — (7, Z) — h + X, (1) 

which describes the inclusive production of a single charged hadron, h. Here, h may 
either refer to a specific charged-hadron species, such as n^, K , or p/p, or to the sum 
of all charged hadrons. The partonic cross sections pertinent to process (|XJ) can entirely 
be calculated in perturbative QCD with no additional input, except for a s . They are 
known at next-to-leading order (NLO) || and even at next-to-next-to-leading order ||. 
The subsequent transition of the partons into hadrons takes place at an energy scale of 
the order of 1 GeV and can, therefore, not be treated in perturbation theory. Instead, 
the hadronization of the partons is described by FF's D^(x, Q 2 ). Their values correspond 
to the probability that the parton a, which is produced at short distance, of order 1/Q, 
fragments into the hadron h carrying the fraction x of the momentum of a. In the case 
of process (|l|), Q is typically of the order of the center-of-mass (CM) energy ^/s. Given 
their x dependence at some scale Q , the evolution of the FF's with Q may be computed 
perturbatively from the timelike Altarelli-Parisi equations , which are presently known 
through NLO TQ. This method to determine af'(Mz) is particularly clean in the sense 
that, unlike other methods, it is not plagued by uncertainties associated with hadroniza- 
tion corrections, jet algorithms, parton density functions (PDF's), etc. We recall that, 
similarly to the scaling violations in the PDF's, perturbative QCD only predicts the Q 2 
dependence of the FF's. Therefore, measurements at different CM energies are needed in 
order to extract values of (M^). Furthermore, since the Q 2 evolution mixes the quark 
and gluon FF's, it is essential to determine all FF's individually. 

In 1994/95, two of us, together with Binnewies, extracted it^ and K FF's through fits 
to PEP and partially preliminary LEP1 data and thus determined af\Mz) to be 0.118 
(0.122) at NLO (LO) J8| (BKK). However, these analyses suffered from the lack of specific 
data on the fragmentation of tagged quarks and gluons to ir^, K^, and p/p hadrons. This 
drawback has been cured in 1998 by the advent of a wealth of new data from the LEP1 and 
SLC experiments [§,|l(],[n|,[l^,[n| . The data partly come as light-, o, and 6-quark-enriched 



samples with identified final-state hadrons (ir ± ,K ± , and p/p) P Jl0| , p4]| or as gluon-tagged 
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three-jet samples without hadron identification fT2| , |T3"lj . This new situation motivates us 
to update, refine, and extend the BKK analysis || by generating new LO and NLO sets 
of 7r ± , K ± , and p/p FF's. By also including in our fits 7r ± , K , and p/p data (without 
flavor separation) from PEP |TJ], with CM energy y/s = 29 GeV, we obtain a handle 
on the scaling violations in the FF's, which allows us to extract LO and NLO values of 
af\Mz)- The latter data [TJ| combines small statistical errors with fine binning in x and 
is more constraining than other data from the pre-LEPl/SLC era. 

The NLO formalism for extracting FF's from e + e~ data was comprehensively described 
in Ref. ||. We work in the MS renormalization and factorization scheme and choose the 
renormalization scale \x and the factorization scale Mf to be p, — Mf = £y/s, except 
for gluon-tagged three-jet events, where we put p = Mf = 2^E jet) with E )et being the 
gluon jet energy in the CM frame. Here, the dimensionless parameter £ is introduced 
to determine the theoretical uncertainty in af'(Mz) from scale variations. As usual, we 
allow for variations between £ = 1/2 and 2 around the default value 1. For the actual 
fitting procedure, we use x bins in the interval 0.1 < x < 1 and integrate the theoretical 
functions over the bin widths as is done in the experimental analyses. The restriction 
at small x is introduced to exclude events in the nonperturbative region, where mass 
effects and nonperturbative intrinsic-transverse-momentum effects are important and the 
underlying formalism is insufficient. We parameterize the x dependence of the FF's at 
the starting scale Q as 

D h a {x,Ql) = Nx a {l-xf. (2) 

We treat N, a, and 8 as independent fit parameters. In addition, we take the asymptotic 
scale parameter A^, appropriate for five quark flavors, as a free parameter. Thus, we 
have a total of 46 independent fit parameters. The quality of the fit is measured in 
terms of the \ 2 value per degree of freedom, Xdfi f° r a ^ selected data points. Using a 



multidimensional minimization algorithm |fLq| , we search this 46-dimensional parameter 
space for the point at which the deviation of the theoretical prediction from the data 
becomes minimal. 

The Xbf values achieved for the various data sets used in our LO and NLO fits may be 
seen from Table [I]. Most of the Xdf values lie around unity or below, indicating that the 
fitted FF's describe all data sets within their respective errors. In general, the Xdf values 
come out slightly in favor for the DELPHI |T0| data. The overall goodness of the NLO 
(LO) fit is given by Xbf = 0-98 (0.97). The goodness of our fit may also be judged from 
Figs. [TJ and |2|, where our LO and NLO fit results are compared with the ALEPH PJT2]] , 
DELPHI 0, OPAL 0, and SLD flj data. In Fig. |, we study the differential cross 
section (l/a to t)da h jdx of process (|l|) for 7r ± , K ± , p/p, and unidentified charged hadrons 
at ^/s = 91.2 GeV, normalized to the total hadronic cross section cr to t, as a function of 



the scaled momentum x = Ip^j^fs. As in Refs. [10,11], we assume that the sum of the 
7r^, K^, and p/p data exhaust the full charged-hadron data. We observe that, in all 
cases, the various data are mutually consistent with each other and are nicely described 
by the LO and NLO fits, which is also reflected in the relatively small Xdf values given in 
Table [l|. The LO and NLO fits are almost indistinguishable in those regions of x, where 
the data have small errors. At large x, where the statistical errors are large, the LO 
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Table 1: CM energies, types of data, and Xdf values obtained at LO and NLO for the 
various data samples. 



s [GeV] 



Data type 



Xdf 



in NLO (LO) 

H 



29.0 



(7 



(J 



<J' 



(all) 
(all) 
(all) 



0.64 
1.86 
0.79 



0.71 
1.40 
0.70 



E 
141 



91.2 



a 



a 



(all) 

(uds) 
(b) ' 



1.28 
0.20 
0.43 



1.40 
0.20 
0.41 



m 
[0 



1.32 (1.44) 



<7 W (all) 

a 77 (uds) 
&" (c) 
&" (b) 



1.28 
0.58 
0.72 

0.57 



1.65 
0.60 
0.73 

0.58 



3.09 
[0 1-87 
1.36 
P| 1.00 



Iff 

[0 0.44 

|0| 0.65 
2.11 

m i.2i 



3.13 
2.17 
1.16 
0.99 



fill 






a 



(all) 



(uds) 
a K (c) 
a* (b) 



0.30 
0.86 
0.53 

0.14 



0.32 
0.79 
0.60 

0.14 



0.45 
0.64 
1.90 
1.23 







in 



a* 



(all) 



£jct [GeV] 



26.2 
40.1 



a p (uds) 
o*(c) 
gP (b) 



0.93 
0.09 
0.11 

0.56 

1.19 
1.03 



0.80 
0.06 
0.14 

0.62 

LT8 
0.90 



10 

JEL 

W 
pi 



0.79 
1.29 
0.92 
0.97 



0.70 
1.28 
0.89 
0.89 
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and NLO results sometimes moderately deviate from each other. In Fig. we compare 
the ALEPH Jl2| and OPAL [[Tj| measurements of the gluon FF in gluon-tagged charged- 
hadron production, with E- ]et = 26.2 and 40.1 GeV, respectively, with our LO and NLO 
fit results. The data are nicely fitted, with Xbf values of order unity, as may be seen from 
Table [IJ. By the same token, this implies that the data [ |T2] , |TB| ] are mutually consistent^. 
The purpose of this letter is to update and improve the determinations of and 



f ] (M z ) from the scaling violations in the FF's. We obtain A§, = 88+^ + 2 3 , MeV at 

(5) 91 o+75 +22 
MS -73-29 



LO and = 2131^ I~ M eV at NLO, where the first errors are experimental and the 

(5) 

second ones are theoretical. The experimental errors are determined by varying Aj^ in 
1 The new FF sets can be obtained from [http://www.dcsy.de/~pocttcr/kkp.htiiLl 
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such a way that the total Xdf value is increased by one unit if all the other fit parameters 
are kept fixed, while the theoretical errors are obtained by repeating the LO and NLO 
fits for the scale choices £ = 1/2 and 2. From the LO and NLO formulas for ai n/ ^(/i) @, 
we thus obtain 

«?>(AW = 0.1181 +0-0068 +0.0006 

3 v ; -0.0069-0.0049 v ; 

af{M z) =0.1170 +0-0055+0.0017 

3 v ; -0.0069-0.0025 ^ h y ! 

respectively. Adding the maximum experimental and theoretical deviations from the 
central values in quadrature, we find A^ = (88±41) MeV and af\M z ) = 0.1181±0.0085 

at LO and A|| = (213 ± 79) MeV and af\hl z ) = 0.1170 ± 0.0073 at NLO. We observe 
that our LO and NLO values of ay\Mz) are quite consistent with each other, which 
indicates that our analysis is perturbatively stable. The fact that the respective values of 
A^ significantly differ is a well-known feature of the MS definition of a,s f \fj) @. 

Our values of and af\M z ) perfectly agree with those presently quoted by the 

Particle Data Group (PDG) Q as world averages, AgjL = 212+^ MeV and af\M z ) = 
0.1185 ± 0.002, respectively. Notice that, in contrast to our LO and NLO analyses, 
the PDG evaluates AgjL from af\M z ) using the three-loop relationship f2|. The PDG 
combines twelve different kinds of af\M z ) measurements, including one from the scaling 
violations in the FF's fl~6| , by minimizing the total x 2 value. The world average cited above 
is then estimated from the outcome by allowing for correlations between certain systematic 
errors. It is interesting to investigate how the world average of af\M z ) is affected by 
our analysis. To this end, we first combine the twelve a^(M z ) measurements reported 
in Ref. |] to find af\M z ) = 0.1181 ± 0.0014 with X 2 = 3.74.0 If we replace the value 
a^'(Mz) = 0.125 ± 0.005 ± 0.008 resulting from previous FF analyses |16|], which enters 



the PDG average, with our new NLO value, then we obtain af\M z ) = 0.1180 ± 0.0014 
with x 2 — 3.21, i.e., the face value of the world average essentially goes unchanged, while 
the overall agreement is appreciably improved. This is also evident from the comparison 
of Fig. 0, which summarizes our updated world average, with the corresponding Fig. 9.1 
in Ref. Jj]]. We observe that the central value of our new NLO result for af\M z ) falls 
into the shaded band, which indicates the error of the world average, while in Fig. 9.1 of 



Ref. M the corresponding central value [16| exceeds the world average by 3.3 standard 



deviations of the latter, which is more than for all other eleven processes. Furthermore, 
our new NLO result has a somewhat smaller error (0.0073) than the corresponding result 
HIE! used by the PDG (0.009). If we take the point of view that our new NLO value of 
a)p(Mz) should rather be combined with the result from the previous FF analyses [JC 



before taking the world average, then the latter turns out to be a ( f\M z ) = 0.1181±0.0014 
with x 2 = 3.29. 

In summary, we presented an updated and improved determination of a^(M z ) from 
the LO and NLO analyses of inclusive light-hadron production in e + e~ annihilation. Our 

2 This result slightly differs from the corresponding one found in Ref. M. 
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strategy was to only include in our fits high-precision LEP1 and SLC data with both flavor 
separation and hadron identification (namely, light-, o, and 6-quark-enri died samples of 
ti ± , K , and pp data) |Tl"|jn| , gluon-tagged three-jet samples with a fixed gluon-jet energy 
12 , 13 1 , and the n ± , K , and pjp data sets from the pre-LEPl/SLC era with the highest 



statistics and the finest binning in x |14]]. Our LO and NLO results for af\Mz) are 
given in Eq. (|3]). They should be compared with the result from scaling violations in 
FF's quoted in Ref. Q, 0.125 ± 0.005 ± 0.008 0]. If we repeat the global analysis of 



Ref. with this result replaced by our new NLO value, then the world average (before 
taking into account estimated correlations between systematic errors) is changed from 
af\M z ) = 0.1181 ±0.0014 with X 2 = 3.74 to af\M z ) = 0.1180 ±0.0014 with X 2 = 3.21, 
i.e., the overall agreement is appreciably improved. 

The II. Institut fur Theoretische Physik is supported by the Bundesministerium fur 
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Figure 1: Normalized differential cross section of inclusive hadron production at y/s = 
91.2 GeV as a function of x. The LO (dashed lines) and NLO (solid lines) fit results are 
compared with data from ALEPH || (triangles), DELPHI |K| (circles), and SLD [11 



(squares). The upmost, second, third, and lowest curves refer to charged hadrons, n^, 
K ± , and p/p, respectively. Each pair of curves is rescaled relative to the nearest upper 
one by a factor of 1/5. 
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x 

Figure 2: Gluon FF for charged-hadron production at Q — 52.4 and 80.2 GeV as a 
function of x. The LO (dashed lines) and NLO (solid lines) predictions are compared 
with three-jet data from ALEPH |~2| with E- set = 26.2 GeV (upper curves) and from 
OPAL [|BJ with Ej Ct = 40.1 GeV (lower curves). The OPAL data and the pertinent 
predictions are rescaled by a factor of 1/100. 
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Figure 3: Summary of the values of af\M z ) from various processes. The errors shown 
represent the total errors including theoretical uncertainties. 
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